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We present phase shift measurements for neutron matter waves in vacuum and in low pressure 
Helium using a method originally developed for neutron scattering length measurements in neutron 
interferometry. We search for phase shifts associated with a coupling to scalar fields. We set 
stringent limits for a scalar chameleon field, a prominent quintessence dark energy candidate. We 
find that the coupling constant /3 is less than 1.9 xlO 7 for n = 1 at 95% confidence level, where n 
is an input parameter of the self interaction of the chameleon held ip inversely proportional to ip n . 

PACS numbers: 95.36.+x 03.75.Dg 


I. INTRODUCTION 

The accelerating expansion of the universe suggests that most of the energy in the universe is ’dark energy’. The 
nature and origin of this energy remain unknown. Candidates for dark energy are either Einstein’s cosmological 
constant or dynamical dark energy, i.e. the so-called quintessence canonical scalar field ip, responsible for the late¬ 
time acceleration of the universe expansion. Chameleon fields are a prime example of dynamical dark energy. Their 
effective mass depends on the energy density of matter in which it is immersed [1]. As a result, in a sufficiently 
dense environment the chameleon field is very massive and, correspondingly, substantially Yukawa-suppressed, i.e. 
very short-ranged. In turn, it is essentially massless on cosmological scales [2, 3]. Because of its sensitivity on the 
environment, such a mass-changing scalar field has been called chameleon. Moreover, the chameleon field always 
couples to matter and generates a fifth force with an effective range inversely proportional to its effective mass. 

All models of dark energy involve a light scalar field 0,1 whose effects on solar system tests of gravity needs to be 
shielded. Three main screening mechanisms [3] have been unraveled so far. The K mouflage and Vainshtein screenings 
are very powerful inside a large domain surrounding the earth, rendering their test in laboratory experiments extremely 
arduous. On the other hand, the chameleon mechanism is at work in the presence of dense objects and can be tested 
in near-vacuum experiments 0 ] . This is the case for the Eotwash 0 and Casimir experiments [6|, where the boundary 
plates are screened. Another way of testing the chameleon mechanism involves small and unscreened objects, like 
neutrons under certain conditions 0. 

Concerning chameleon models, a chameleon-photon coupling g e s = /3 7 /Mpi has been proposed, and the detailed 
analysis of the chameleon -photon interaction and a comparison with the cosmological data has been carried out in 0- 
H~2| . A search for photon-chameleon -photon transition has been performed by the experiment CHASE (the GammeV 
CHameleon Afterglow SEarch)[l3| and by the Axion Dark Matter experiment (ADMX) flUl. A search for chameleon 
particles created via photon-chameleon oscillations within a magnetic field is described in [l5|. 

Searches with neutrons directly test the chameleon-matter interaction [3 and do not rely on the existence of a 
chameleon-photon-interaction. The coupling /3 is restricted from below, e.g. (3 must be larger than 50 at n = 1 jl6| . 
and experiments with neutrons have the potential ultimately to find a chameleon field or exclude it in the whole 
parameter space. 

As it has been pointed out by Pokotilovski na , the use of a neutron Lloyd’s interferometer for measurements of the 
phase-shift of the wave function of cold neutrons should allow to determine the chameleon-matter coupling constant. 
The ^Bounce collaboration has searched for the chameleon field using gravity resonance spectroscopy and ultra-cold 
neutrons [l8l - (2l| . In a recent experiment [22j], the upper limit for (3 has been determined as f3 < 5.8 x 10 8 which is 
five orders of magnitude below the previous limit determined by atomic spectra h(> . 

Here we present a new search for chameleon fields by means of neutron interferometry as proposed in 0. The 
self interaction of the chameleon field ip and its interaction to an environment with mass density p are described by 
the effective potential [2^, l24| 
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where /3 is the coupling constant, n is an input parameter (the so-called Ratra-Peebles index) and A ss 2.4 x 10 12 
GeV defines an energy scale Q. Mpi = ^Hc/(8ttG) = 4.341 x 10~ 9 kg denotes the reduced Planck mass. The 
chameleon field ip creates a potential for neutrons given by V = /3 <pm/Mp\ where m denotes the neutron mass. When 
passing this potential, neutrons accumulate the phase 

c = S v{x)ix = ~w j B lk l ' e{ - x)dx ' <2) 

where k denotes the neutron wave vector modulus k = 27r/A. 

For strong coupling (/? 1) the chameleon field is suppressed at the presence of matter, even at low mass densities 

like air at ambient pressure. Only in vacuum the chameleon field can persist. By placing a vacuum cell into one arm 
of the neutron interferometer and allowing ambient air in the other arm we can directly probe the chameleon field. 
The setup resembles a standard setup for measuring neutron scattering lengths [HI, but instead of measuring the 
phase shift of sample material we measure the phase shift of vacuum. 

The chameleon field vanishes at the walls of the vacuum chamber but increases bubble-like towards the middle of 
the chamber, cf. Fig. |T](c). The more the field increases the lower the remaining gas pressure is, i.e. the better the 
vacuum is. Thus we have two options of performing a relative phase measurement which is necessary to cancel the 
unknown intrinsic interferometer phase and the air phase shift. In the pressure mode we vary the pressure in the 
vacuum cell by letting in different amounts of Helium. In the profile mode we keep the pressure constant but move 
the chamber transversally to the beam in order to record a profile of the chameleon bubble. Neither method detects 
any chameleon-like signature, giving rise to new constraints of the chameleon theory. 


II. SETUP 

The experiment is carried out at the neutron interferometry setup S18 at the Institut Laue-Langevin (ILL) in 
Grenoble. A perfect crystal silicon interferometer is used, Fig. Q](a), at 45° Bragg angle and 2.72 A mean wave length 
A with 0.043 A wavelength distribution width (FWHM). The two beam paths within the interferometer are separated 
by 50 mm over a length of 160 mm. Neutron detectors with an efficiency above 99% measure the intensities of the two 
exit beams labeled O and H respectively. A vacuum chamber with inner dimensions 40 x 40 x 94 mm is inserted in 
the left or right beam path. The other beam path always contains one of the two air chambers which sit alongside the 
vacuum chamber. The whole chamber box can be moved sidewards for swapping the vacuum cell between the left and 
the right beam path and to probe different beam trajectories within the vacuum cell. The air chambers ensure that 
both beam paths contain the same amount of wall material (aluminium). In addition, the extension of the vacuum 
cell by air chambers minimizes possible disturbances of the thermal environment of the crystal when the chamber box 
is moved. We label different chamber positions by the letters ’a’ to ’n’ as indicated in the figure. 

The air chambers are connected to ambient air by a hole in the top of the chambers. The vacuum chamber is 
connected to a vacuum control system consisting of pressure gauge, motorized leak valve and pumps, as indicated in 
Fig. CO (b). The pumps (pre-pump and turbomolecular pump) are running continuously while a controlled amount of 
Helium is let in through the leak valve in order to tune the pressure. The pressure gauge is corrected for the use with 
Helium. 


III. DATA ACQUISITION AND EVALUATION 

Phases in neutron interferometry are measured by rotating an auxiliary phase flag and recording the intensity 
oscillations between O and H detector, cf. Fig. [0 Such interferograms are measured before and after some parameter 
change. The shift of the sine curves with respect to each other represents the phase shift induced by the parameter 
change. The recording of each interferogram takes typically half an hour, and during that time the intrinsic phase of 
the interferometer can drift due to temperature changes or other environmental factors. To compensate such drifts 
we interlace phase flag movement and parameter change. The phase flag is rotated to the first angular position and 
neutrons are counted for a certain amount of time for each parameter setting. Then the phase flag is rotated to the 
next position and neutrons are counted again for all parameter settings etc. In the end we obtain interferograms 
measured simultaneously for all parameter settings. Their relative phases are free of phase drifts. 

We use the largest neutron interferometer available [26| with a loop size of 50 x 160 mm in order to maximize the 
size of the vacuum cell. Such big single crystal interferometers are extremely sensitive to temperature gradients, air 
flow, vibrations, bending, etc. Hence the interference contrast (fringe visibility) is restricted to about 10% to 30%. 
The interferograms look a bit more noisy than what can be explained by pure counting statistics. This means that 
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FIG. 1: (a) Top view of the interferometry setup shown in chamber position ’h’. The chamber box (blue) can be moved 
transversally allowing the beams to pass at different positions, labeled by ’a’ to ’n’. (b) Scheme of the vacuum handling and 
axial view of the vacuum chamber, (c) Longitudinal and transverse bubble shape of the chameleon field in the vacuum cell. 
The beam positions ’a’, ’d’ and ’g’ are indicated by green rectangles. 



FIG. 2: Recorded intensity oscillations between O and H detector as a function of the optical path length difference AD created 
by rotating the phase flag. The three curves in red, purple and blue represent the interferograms at the ’a’, ’d’ and ’g’ position 
respectively. The phase shift between these raw curves is created by position dependent wall thickness variations, cf. Fig. [3] 


the phase is slightly fluctuating within the recording time of each interferogram. We conservatively account for this 
noise by performing a y 2 test for each sine fit and by blowing up the fit error (by a factor of about 2) such that the 
y 2 test is satisfied. 
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FIG. 3: (a) Mechanical measurement: measured total thickness of the entry and exit walls at two different vertical positions. 
The thickness is increased around the screw holes (black dots). The positions of beam path I and II are indicated for several 
chamber positions, (b) Phase shift calculated from (a), caused by the different amount of wall material (aluminium) in path 
I and II for different chamber positions, (c) Phase profile of the vacuum chamber with and without correction for the wall 
thickness. On the left side (a-g) path I passes air and path II passes vacuum; vice versa on the right side. 


A. Profile mode 

In the profile mode we measure at up to 14 transverse beam positions for each phase flag position in order to 
look for bubble-like phase profiles. Ideally, the entry and exit walls of the vacuum chamber are flat and parallel and 
therefore would not alter the phase if the chamber is transversally moved. Unfortunately, the screw holes of our walls 
have been drilled after the surfaces had been polished. As a consequence, the surfaces are elevated by a few microns 
around each screw hole, and all our data in profile mode require a position dependent phase correction, based on a 
careful mapping of the wall thickness, as indicated in Fig. [3] (a) and (b). The screw hole positions are indicated by 
black dots in Fig. [3] (a) and Fig. |T](b). Fortunately, the correction depends mainly linear on the beam position, while 
the shape of the chameleon bubble is expected to be mainly parabolic. To be precise, the chameleon profile must be 
symmetric with respect to the cell center, and can therefore consist only of quadratic and higher even orders. Thus, 
there is no danger that the wall thickness correction completely mimics or hides the chameleon feature. We determine 
the thickness correction at the upper and the lower edge of the beam, and use the average as correction and a quarter 
of the difference as uncertainty of the correction. 

Fig. [3] (c) shows the recorded phase over a complete profile. In position ’a’ to ’g’ path I passes air and path II 
passes vacuum while in position ’h’ to V it is the other way round, cf. Fig. |T] (a). The slope within each group is 
caused by the thickness variation of the chamber walls. The step between the two groups comes from the sign change 
of the air phase shift when air and vacuum are swapped between the beam paths. 

The height of the chameleon bubble can be determined by comparing the phase at the center of the chamber with 
the phase at its side, close to the chamber walls. Therefore we make most of our measurements at positions ’a’, 
’d’ and ’g’. Fig. []] (a) summarizes the result of the bubble height measurements for various pressure settings. The 
statistical error of the phase can be reduced to typically 0.9° by averaging over 15 measurements. However, the 
thickness correction, which is applied after the statistical averaging, increases the error again to typically 2.5°. 
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(a) phase at center minus phase at side 



(b) phase at low pressure minus phase at 0.011 mbar 
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FIG. 4: Measured phase shifts in the profile mode (a) and pressure mode (b) compared to calculations for different values of 
P and n. (c) Exclusion plot comparing our results with other experiments. The limit for p at 95% confidence level is shown 
for different values of n. 
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B. Pressure mode 

In the pressure mode we apply four different pressures at each phase flag position. A quick pressure change is only 
possible in the pressure range of the turbomolecular pump, i.e. below 10“ 2 mbar. The average of four such runs is 
shown in Fig. [I](b). In order to compensate phase drifts we use the phase at the highest pressure (0.011 mbar) as 
reference and determine the phase shift between this pressure and the other pressure values. The magnitude of the 
phase shift created by the Helium itself is in the order of < 0.001° in this pressure range and can be neglected. 


IV. LIMIT CALCULATION 


The solution of the chameleon field in vacuum confined between two walls at x = 


±d/2 is given analytically [13 by 


<Pm(x) = A 
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For higher dimensions and for finite gas pressure the calculation has to be done numerically. The left side of Fig. |T| 
(c) shows the longitudinal field profile <pzo{x, y = 0, z = 0) along the center of the chamber calculated in 3D and for 
vacuum. It vanishes at the walls and increases towards the middle. Over most of the range it is nearly constant because 
it is limited by the much narrower transverse confinement. The transverse field distribution (p 2 D{y,z) ~ <^ 3 d(0, y, z) 
is shown on the right side. 

Since the full 3D calculation is very time consuming we assume to good approximation that the field depends only 
on the transverse coordinates. We account for the longitudinal drop close to the walls by calculating an effective 
chamber length 4ff such that 


A/2 


1 - 1/2 
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Thus the true length l = 94 mm reduces effectively to Z e ff = {84,85.6,86.8,87.6} mm respectively for n = {1, 2, 3,4}. 
The expected phase shift £ given by Eq. 0 simplifies to 
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and is plotted in Fig. [I] (a) and (b) assuming various values for /3 and n. 

We calculate limits for /? by comparing the calculated phase shifts £ with the measured phase shifts ( ± a. We 
assume certain values of /? and calculate the corresponding probability p. 
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The sum goes over the data points shown in Fig. 0] (a) and (b). We determine the limit Aim with 95% confidence 
level by numerically solving the equation 



p{P) d/3 = 95%. 


The calculation is repeated for different values of n yielding the following results. 


Aii 


1.9 x 10 7 , n = l, 

5.8 x 10 7 , n = 2, 
2.0 x10 s , 7i = 3, 

4.8 x10 s , 7i = 4. 


( 8 ) 


( 9 ) 


Both the profile mode and the pressure mode contribute about equally to these limits. For large n the profile mode 
becomes less sensitive because the bubble shape becomes flatter on the top. 


V. CONCLUSION 

Our search for chameleons by means of neutron interferometry failed in finding ones but succeeded in deriving new 
upper bounds for the coupling constant /3, listed in Eq. 0. For n = 1 the new limit is a factor of 30 below the 
previous one which has been obtained by gravity resonance spectroscopy, cf. Fig. 01(c). There remains a range of five 
orders of magnitude for j3 where chameleons have not been excluded yet. 
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Note added - During the preparation of this Letter, another interfometric experiment searching for the chameleon, 
using ultracold Cs atoms, has been reported f2S 
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